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A.  Problem  Statement 

This  project  concerned  numerical  methods  for  solving  linear  systems  of  equations  of  the  type 
arising  from  discretization  of  boundary  value  problems  of  elliptic  and  mixed  type.  The  problems 
considered  were  of  fundamental  use  in  mathematical  models  used  in  structural  analysis  and 
fluid  dynamics.  The  emphasis  was  on  preconditioning  techniques  in  which  properties  .of  the 
problem  are  used  to  construct  approximations  that  are  easy  to  compute  with  and  that  lead  to 
rapid  convergence  of  iterative  methods.  The  approach  included  both  analytic  studies  producing 
bounds  on  convergence  rates,  and  computational  experiments  that  confirm  and  supplement  the 
analysis. 

B.  Summary  of  Major  Results 

1.  Iterative  methods  for  the  discrete  convection-diffusion  equation.  An  analysis  has 
been  conducted  of  the  effect  of  ordering  of  unknowns  on  performance  of  relaxation  methods 
for  solving  linear  systems  arising  from  discretization  of  the  (steady-state)  convection- diffusion 
equation.  Main  results  consist  of  analytic  bounds  on  norms  of  relaxation  operators  associated 
with  both  one-dimensional  and  two-dimensional  problems.  The  results  rigorously  establish  that 
relaxation  that  “follows  the  flow”  associated  with  the  problem  displays  fast  convergence,  but 
that  there  are  latencies  associated  with  relaxation  based  on  red-black  ordering  and  sweeping  in 
a  direction  opposite  to  the  direction  of  flow. 

2.  Parallel  algorithms  for  discrete  elliptic  problems.  A  study  of  parallel  implementation  of 
the  /ip- version  of  the  finite  element  method  on  a  shared-memory  computer  shows  that  performance 
can  be  improved  by  grouping  computations  into  blocks  to  avoid  memory  conflicts.  This  enhances 
performance  on  shared  memory  machines  and,  more  importantly,  is  needed  to  make  efficient  use 
of  hierarchical  memory  on  arbitrary  architectures. 

3.  Solution  algorithms  for  the  Stokes  and  Navier-Stokes  equations.  The  most  recent 
emphasis  has  been  on  numerical  solution  techniques  for  problems  arising  in  models  of  incom¬ 
pressible  fluid  flow.  The  typical  linear  systems  produced  are  nonsymmetric  and  indefinite,  which 
present  challenges  for  solution  algorithms,  but  the  structure  of  the  problems  can  be  used  to 
develop  effective  algorithms. 

a.  New  preconditioning  techniques  for  the  linearized  discrete  Navier-Stokes  equations  have 
been  developed  that  produce  preconditioned  linear  systems  with  eigenvalues  that  are  bound¬ 
ed  independent  of  the  discretization  mesh  size.  Experimental  results  indicate  that  Krylov 
subspace  iterative  methods  combined  with  these  preconditioners  display  rates  of  conver¬ 
gence  that  are  independent  of  the  mesh  parameter.  In  addition,  they  show  that  performance 
can  be  improved  through  the  use  of  inner  iterations  for  solving  subproblems  consisting  of  a 
set  of  discrete  convection-diffusion  equations. 

b.  An  analysis  and  empirical  study  of  a  variant  of  the  Uzawa  method  for  self-adjoint  saddle 
point  problems  shows  that  the  most  costly  part  of  the  algorithm,  solution  of  a  symmetric 
positive- definite  subsystem  at  every  iteration,  can  be  replaced  by  an  approximate  solution 
obtained  by  an  inner  iteration.  The  convergence  behavior  of  the  resulting  “inexact  Uzawa” 
algorithm  is  relatively  insensitive  to  the  accuracy  of  the  approximate  solution.  Analysis 


2 


and  experiments  with  the  Stokes  equations  confirm  and  supplement  these  results. 

c.  An  empirical  comparison  of  four  solvers  for  the  discrete  Stokes  equations  (variants  of  the 
Uzawa,  preconditioned  conjugate  gradient,  preconditioned  conjugate  residual,  and  multi¬ 
grid  methods)  indicates  that  multigrid  with  smoothing  based  on  incomplete  factorizaton 
is  more  efficient  than  the  other  methods,  but  typically  by  no  more  than  a  factor  of  two. 
The  conjugate  residual  method  has  the  advantages  of  being  both  independent  of  iteration 
parameters  and  widely  applicable. 

4.  Other  results  for  elliptic  problems.  A  comparison  of  the  performance  of  low  order  and 
high  order  finite  elements  for  solving  the  discrete  Poisson  equation  has  been  made,  with  results 
showing  the  efficiences  achieved  by  high  order  hp-elements. 
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